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Abstract 



Amorphous materials as diverse as foams, emulsions, colloidal suspensions and granular me- 
dia can jam into a rigid, disordered state where they withstand finite shear stresses before yield- 
ing. The jamming transition has been studied extensively, in particular in computer simulations 
of frictionless, soft, purely repulsive spheres. Foams and emulsions are the closest realizations 
of this model, and in foams, the (un)jamming point corresponds to the wet limit, where the 
bubbles become spherical and just form contacts. Here we sketch the relevance of the jamming 
perspective for the geometry and flow of foams — and also discuss the impact that foams studies 
may have on theoretical studies on jamming. 

We first briefly review insights into the crucial role of disorder in these systems, culminating 
in the breakdown of the affine assumption that underlies the rich mechanics near jamming. 
Second, we discuss how crucial theoretical predictions, such as the square root scaling of contact 
number with packing fraction, and the nontrivial role of disorder and fluctuations for flow have 
been observed in experiments on 2D foams. Third, we discuss a scaling model for the rheology 
of disordered media that appears to capture the key features of the flow of foams, emulsions 
and soft colloidal suspensions. Finally, we discuss how best to confront predictions of this model 
with experimental data. 



1 Introduction 



Amorphous materials as diverse as foams, emulsions, colloidal suspensions and granular media can 



jam into a rigid, disordered state where they withstand finite shear stresses before yielding [Liu and 
Nagel I p^98| ); |0'Hern et al. | pOOS] ); |Liu and Nagel | ( [Mol ); [van Hecte] ( [Mol )]. The jamming 
transition has been studied extensively, in particular in computer simulations of frictionless, soft, 
purely repulsive spheres. Foams and emulsions are the closest realizations of this model, and in 
foams, the (un)jamming point corresponds to the wet limit, where the bubbles become spherical 
and just touch. What does jamming tell us about foams? And what do experiments on foams teach 
us about jamming? 

Liquid foams and emulsions are dispersions of gas bubbles or droplets in a immiscible (sec- 
ond) liquid phase, stabilized by surfactants. Phenomenologically, these materials exhibit plastic 
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flow under large stresses, but behave elastically under small stress - their macroscopic response to 
mechanical perturbations is a complex mix of elastic, plastic and viscous efl'ects fHohler and Cohen- 



Addad (2005)] - as is the case for many other disordered soft materials ranging from colloidal 



suspensions to granular materials Larson (1999)]. 

How should we think about the mechanics of such materials, and more specifically, about the 
mechanics of foams? More precisely, how are the laws that govern the local dynamics of bubbles and 
soap films related to the macroscopic, collective behavior of a foam? It is useful to introduce two 
extreme points of view here first, being fully aware that most researchers' views are more subtle. 

One extreme point of view would be that a complete knowledge of the local mechanics is sufficient 
to capture the global behavior - from this perspective, the translation from local to global behavior is 
nothing but a straightforward coarse-graining procedure, and nontrivial collective behavior, such as 
nonlinear rheology, always finds its root cause in similar nontrivial local events, such as a nonlinear 
viscous friction law. For reasons that will become clear, we refer to this as the 'affine' picture. 

Another extreme point of view would be that none of the local details matter and that the global 
behavior of disordered media, and foams in particular, is set by "universal" collective mechanisms 
[Gates et al.~| p997l ); psson and TeTteT] ^So7\) ; [Argon | ( |T979| ); |Falk and Langer | ( |T998| )]. Hence 
knowledge of the detailed interactions is irrelevant, and different local interactions, such as different 
local friction laws, may lead to the same macroscopic behavior, i.e. rheology. 

Here we put forward a third point of view: qualitatively new behavior can emerge at the global 
level, but local details still may matter. So, in our view, there is no universality in the sense that 
changes in the microscopies often lead to changes in the macroscopic properties. However, due 
to the general 'non-affine' and strongly fiuctuating nature of the mechanics of disordered systems, 
the translation from the local to the global level can be also highly nontrivial. We suggest that 
the mechanisms that connect the micro and macro world are universal. As an example of such a 
robust mechanism, we will discuss the role of the strength of the dynamic fiuctuations, which in 
turn depend on the wetness of the foam and the fiow rate. 

We will illustrate our point of view by discussing recent experimental and numerical work on 
foams, and by making explicit contact with the jamming framework that has been emerging in 
recent years. Jamming refers to the creation (or loss) of rigidity in disordered systems in general, 
but the most studied models in jamming are close to models for wet foams and actually were inspired 
by earlier foam work |Bolton and Weaire (1990); Durian (1995)]. 

The outline of this paper is as follows. In section 2 we discuss the basic jamming scenario for 
foams, including the differences between the mechanics of disordered foams and ordered models 
— we conclude that disorder leads to strong non-affinity, which becomes increasingly dominant 
near jamming. In section 3 we discuss our recent observations of the nontrivial scaling law that 
relates contact number and packing fraction of foams and jammed materials — this scaling law was 
first seen in simulations twenty years ago Bolton and Weaire (1990)], and has now finally been 
observed in 2D foams [Katgert and van Hecke (2010)]. Section 4 deals with our recent observations 
of the difference in rheology between ordered and disordered foams, and evidences both a nontrivial 
rheology for disordered 2D foams and a nontrivial scaling of the dynamic fiuctuations |Katgert~et 



al. (2008); Mobius et al. (2010)]. In section 5 we sketch the outlines of a model for the rheology of 



disordered media near jamming that we recently introduced Tighe et al. (2010)] which reproduces 
the experimental observations discussed in section 4 well. In section 6 we briefly outline how best to 
confront our predictions with experimental data, and we close with a short discussion and outlook. 
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2 Jamming and Non-affinity: Consequences of Disorder 



Jamming Scenario for Foams — Some of the earliest studies that consider the question of 
rigidity of packings of particles concern the loss of rigidity in foams and emulsions with increasing 

'he gas fraction (j) plays 



(1986 


); 


Kraynik ( 


1988) 



a crucial role in determining the foam's structure and rigidity. The interactions between bubbles 
are repulsive and viscous, and static foams are similar to the frictionless soft spheres used in most 
models for jamming. In real foams, gravity (which causes drainage) and gas diffusion (which causes 
coarsening) may play a role, although these effects can be minimized by studying quasi-2D foams 



and using inert gases Weaire and Hutzler (1999)]. 

The (un)jamming scenario for foams is illustrated in Fig. 1. When the gas fraction approaches 
one, the foam is called dry. Application of deformations causes the liquid films to be stretched, and 
the increase in surface area then provides a restoring force: dry foams are jammed. When the gas 
fraction is lowered and the foam becomes wetter, the gas bubbles become increasingly spherical, 
and the foam loses rigidity for some critical gas fraction (j)c where the bubbles lose contact (Fig. 1). 
The unjamming transition is thus governed by the gas fraction, which typically is seen as a material 
parameter. For emulsions the same scenario arises [Mason et al. (1995)]. 



Vacuum 



Marginal 



Rigid 



Jamming 




Rigidity loss 



Figure 1: Topviews of 2D foams, consisting of a mix of 2 and 3 mm bubbles trapped below a top 
plate. At low packing fractions (left), the bubbles do not form contacts and the materials is in a 
mechanical vacuum state. At high packing fractions (right), the bubbles are squeezed together and 
form a jammed, rigid state. At intermediate packing fractions, the bubbles just touch and form a 
marginal state. 



Elasticity of Disordered Media — Foams and emulsions typically form highly disordered 
packings. Is this disorder important for the mechanical properties, and if so, how can one deal with 
it? 

The simplest approach is to ignore disorder altogether and start from ordered, "crystalline" 
packings. Tacitly assuming that all bubbles experience the same deformation, analytical calculations 
are then feasible, because one only needs to consider a single bubble and its neighbors to capture 
the foams' geometry and mechanical response. 

A famous example are two-dimensional hexagonal packings of monodisperse bubbles ("liquid 
honeycombs" [Princen" (1983); Kraynik (1988)]). The only control parameter is then the liquid 
fraction of the foam or, alternatively, the packing fraction of the gas bubbles </> (the sum of liquid 
fraction and packing fraction is one). 

The most striking features of these liquid honeycombs are as follows, (i) The bubbles lose 
contact at the critical density 0c equal to -^^^ ~ 0.9069, and ordered foam packings are jammed 
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for larger densities, (ii) When for such a model foam (p is lowered towards (j)c, the yield stress and 
shear modulus remain finite and of similar order, and jump to zero precisely at (pc- (Hi) The contact 
number (average number of contacting neighbors per bubble) remains constant at 6 in the jammed 
regime. Similar results can be obtained for three-dimensional ordered foams, where (pc is given by 
the packing density of the HCP lattice ~ 0.7405 |Hohler et al. 



3v^ 



(2008)]. 



Confronted with experiments on foams and emulsions, or numerical experiments on disordered 
foams, all three of these predictions fail! First, the critical packing fraction is substantially lower. 



around 84% in 2D and 64% in 3D Bolton and Weaire ( 1990 ); Durian ( 1995 ); Mason et al. ( 1995 ); 



Saint- Jalmes and Durian (1999); Lespiat et al. (2011)]. The fact that the critical packing density 



for ordered systems is higher than that for disordered systems may not be a surprise, given that 
the bubbles are undeformed spheres at the jamming threshold, and it is well known that ordered 
sphere packings are denser than irregular ones [Weitz (2004)]. 

Second, the yield stress and shear modulus vanish smoothly when the critical packing fraction 
is reached. Early evidence comes from measurements for polydisperse emulsions by [Princen and 
Kiss (1986), who observed a substantial lowering of the shear modulus when cj) is lowered. Later 



measurements by Mason et al. (1995, 1996, 1997) and Saint- Jalmes and Durian (1999) of the 



shear modulus and osmotic pressure of compressed disordered emulsions and foams found similar 
behavior for the loss of rigidity. When scaled appropriately with the Laplace pressure, which sets 
the local "stiffness" of the droplets, the shear modulus grows continuously with cj) and vanishes at 
(j)c ~ 0.635, which corresponds to random close packing in three dimensions. There is also ample 
numerical evidence for this [Bolton and Weaire (1990); Durian (1995); O'Hern et al. (2003); van 
Hecke (2010)] - we will come back to this below. 

Third, the contact number was found to vary smoothly with packing fraction. Early evidence 
comes from numerical simulations by Bolton and Weaire (1990), who numerically probed how a 



disordered foam loses rigidity when its gas fraction is decreased . It was found that the contact 
number z decreases smoothly with (p. At (/> = 1 the contact number equals six, as expected. When 
(f) —7- (j)c, the contact number appears to reach the marginal value, Zc = 4. In related work on the 
so-called bubble model developed for wet foams in 1995, Durian reached similar conclusions for 
two-dimensional model foams, and found that the contact number indeed approaches 4(= 2D) near 
jamming, and observed the nontrivial square root scaling of z — 4 with excess density for the first 
time. All these findings are consistent with what is found in closely related models of frictionless 



soft spheres near jamming O'Hern et al. 



(2003 


); 


van Hecke 


(2010) 



We thus conclude that disorder plays an absolutely crucial role, and in general cannot be treated 
as a perturbation from the regular, ordered case. 

Non-affinity — In the preceding subsection we presented a wealth of simulational and exper- 
imental evidence that invalidates simple predictions for the elasticity of disordered media based on 
intuition derived from ordered packings. The crucial ingredient that is missing is the non-affine 
nature of the deformations of disordered packings (Fig. [2]). In an affine deformation, all the local 
deformations simply follow the global deformation, implying that the local motion of the particles 
is as if they where pinned to a rubber sheet. In a disordered system such as a foam, where particles 
only interact with their neighbors, the local, disordered environment is crucial. The key observation 
is that, while materials far away from the jamming point (Fig. 2a) exhibit deformations that are 
close to affine (Fig. 2b), materials closer to the jamming point (Fig. 2c) exhibit deformations that 
become increasingly disordered and non-affine (Fig. 2d). 

The role of these non-affinities is to critically change the elastic shear modulus of jammed 
materials such as foams O'Hern et al. (2003)]. The scenario is as follows: far away from jamming, 
the local deformations are similar to the global ones - so a 1% shear strain leads to a typical change 
of the local deformations that is also of the order of 1%. Since these changes of deformations are 
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Figure 2: (a) Dense numerical bubble packing, (b) Shear response of packing depicted in (a): 
the bubbles move affinely and thus largely follow the imposed strain, (c) Bubble packing close to 
jamming and its shear response (d): The bubbles exhibit non-affine motion and swirly flow patterns. 



associated with changes in the forces between bubbles/particles, and concomittant changes in the 
elastic energy, one can then immediately estimate the change in elastic energy once the interactions 
between particles are known. In particular, for linear interactions with a local spring constant k, 
one concludes that the shear modulus G on the order of k, just as in the afRne picture sketched 
above. 

Closer to jamming, however, the number of contacts per bubbles drops, and bubbles are increas- 
ingly free to deviate from the afRne field so as to minimize the changes in elastic energy, and, as a 
consequence, the shear modulus G gets smaller than would be expected from an affine assumption. 
At the jamming point, the shear modulus vanishes for harmonically interacting particles — clearly 
here the translation from local to global interaction is nontrivial! For different (e.g., power law in- 
teractions such a Hertzian interactions), the shear modulus behaves different, so there is no simple 
universality. However, the ratio of shear modulus to local spring constant (or equivalently, the ratio 
of shear modulus to bulk modulus K) scales in a robust manner with distance to jamming: indepen- 
dent of dimension and interaction potential, G/K oc z — Zc, where Zc is the number of contacts at the 
critical point (=2D). [O'Hern et al. (2003); Ellenbroek et al. (2009)]. This illustrates our earlier 
point: local interactions matter, but there are universal and sometimes nontrivial mechanisms that 
translate the local interactions to the global behavior. 

As there is no simple way to estimate the particle motions and deformations in disordered 
systems, one needs to resort to (numerical) experiments. Jamming can be seen as the avenue 
that connects the results of such experiments. Jamming aims at capturing the mechanical and 
geometric properties of disordered systems, building on two insights: first, that the non-affine 
character becomes dominant near the jamming transition, and second, that disorder and non- 
affinity are not weak perturbations away from the ordered, affine case, but may lead to completely 



new physics 


Somfai et al. 


( 


2005 


); Mason et al. (1995 


); 


Radjai and Roux 


(2002 


2004 


); 


Lemaitre and Maloney ( 


2006 


); Maloney 


(2006 


)]• 



3 Contacts and density in 2D foams near jamming 



We have recently experimentally investigated the static structure of a disordered monolayer of foam 
bubbles that form a bidisperse packing similar to the detail shown in Fig. 3. This monolayer floats 



on the surface of a soapy solution and is bound on the top by a well-leveled glass plate Katgert 



et al. (2008); Katgert and van Hecke (2010)]. This setup allows for direct optical access of all 
bubbles, in contrast to three dimensional foams, which are opaque due to multiple scattering of the 
interfaces. The packing fraction can be varied simply by in- or decreasing the gap between the glass 
plate and the soapy solution. 

The question we set out to answer is whether the numerically ubiquitous scaling of the excess 



contact number with the square root of the excess density |Durian (1995); O'Hern et al. (2003)] 



could be observed in experiment. In order to probe this question we prepare many distinct packings 
of our bidisperse foam at fixed (j), by stirring the bulk soapy solution to rearrange the packing. After 
the packings have relaxed to an equilibrium state, we take photographs of the resulting structure 
with a 6 megapixel camera and analyse these by advanced image analysis. Even though the bubble 
are three dimensional entities, we adjust the lighting of our experiment such that we image the 
bubbles as two-dimensional discs, at the point where they are broadest. This approach seems to 
be justified, a posteriori. For each realization we determine the average contact number z and the 
packing fraction 0, see grey dots in Fig. 3. For high packing fractions, the contact number tends to 
z = 6, the expected value for disordered cellular structures [Weaire and Hutzler (1999)]. However, 
we observe that, as the foam packing fraction is reduced, the average contact number decreases, 
ultimately reaching Zc = 4, at a packing fraction around (p = 0.84. Below this value the foam loses 
stability. Both the value of the critical density and contact number at the unjamming point are 
fully consistent with predictions. 

We also plot the average over all images for each packing fraction (black circles). Clearly the 
variation of z with (f> is similar to a square root, and to compare our data to predictions we fit 
our data to a power law fit of the form z = 4: + zq * {(p — (f>c)^ (red curve in Fig. 2. The best fit 
gives us zq = 4.02 ± 0.20, (pc = 0.842 ± 0.002 and /3 = 0.50 ± 0.02, in remarkable agreement with 
theoretical predictions by O'Hern et al and Durian [O'Hern et al. (2003); Durian (1995)], who 
found zo = 3.6 ± 0.5, (pc = 0.841 ± 0.002 and /3 = 0.49 ± 0.03. 

In simulation studies the packing fraction is often calculated by counting the overlaps between 
bubbles twice. In the experiment we are restricted to extracting an area fraction from images. We 
therefore convert our experimentally determined packing fractions to the corresponding theoretical 
packing fractions. If we do not convert our experimentally determined packing fractions, we find 
an effective scaling exponent /3 = 0.70 when plotting z as a function of pexp , see inset of Fig. 3. 

We were not the the first to experimentally investigate the scaling of z with (p. Majmudar et 
al. Majmudar et al. (2007)] have extracted the same quantities from images of two-dimensional. 



frictional, photoelastic discs and compared these to predictions from simulations. From their data 
it appears the prefactor Zq ~ 16, inconsistent with simulations. Our results do allow for a direct 
comparison with predictions for frictionless jamming, which can be seen from the excellent agree- 
ment between parameters. We also note that Brujic et al. recently have probed this squareroot 
scaling in 3D emulsions (private communication). 

We conclude that, as far as static structure is concerned, (2D) foams and numerical models for 
jamming are in excellent agreement. Disorder is crucial for real foams. 



6 



N 



6.0 - 




0.84 



0.84 0.88 , 0.92 0.96 
I 

exp 



0.90 



0.96 



1.02 



1.08 



th 



Figure 3: Average contact number versus (j) for experimental bidisperse foams: grey dots indicate 
data for each individual realization, black circles indicate averages for each globally set packing 
fraction. Solid red line is fit of the form z = 4+ zq *{(/)- cl)cf , with zq = 4.02±0.20, (/)c = 0.842±0.002 
and P = 0.50 ± 0.02. Inset is data plotted versus experimentally determined packing fraction c/'exp- 
The fit has a power law exponent of 0.70. 



4 Disorder and fluctuations in foam flows: Experiments 



In a series of experiments |Katgert et aT~| ( [20081 [20091 ) ; [Mobius etaFlpOlOl )] 

we have recently probed 

the role of disorder in the flow of 2D foams. In these experiments, we shear a foam monolayer that 
is floating on a soapy solution and bound by a glass plate. We calculate averaged velocity profiles 
and track the bubble displacements. Our two main findings are 1) that disorder plays a crucial 
role in the rheology — the global rheology is very different than what you would expect based on 
measurements of the local drag forces — and that 2) fluctuations become increasingly strong as the 
strain rate is lowered. 

Drag forces — To probe the role of disorder, we have compared the averaged flow profiles 
and bubble trajectories of a ordered, crystalline foam and a disordered, bidisperse foam (Fig. 4). 
For monodisperse foams (Fig. 4a), the bubbles move past each other in a ziplike fashion along the 
crystal planes of the hexagonal lattice formed by the bubbles, see Fig. 4b for tracks. The velocity 
profiles (Fig. 4b) are strongly localized and are independent of driving velocity. For disordered foams 
(Fig 4c)the situation is vastly different: the bubble motion becomes disordered and organizes into 
the swirly, collective patterns typically seen in materials near jamming, see Fig 4d. The velocity 
profiles now depend on the driving velocity, and become increasingly delocalized as the driving 
velocity is decreased. 

To understand the flow profiles, we introduce a simple force balance model in which the averaged 
viscous friction between bubbles that move past each other, F^, is balanced by the viscous friction 



due to the bubbles moving past the glass plate, Janiaud et al. (2006); Katgert et al. (2008 



2009) 



We note here that without the top plate, the flow profiles are essentially linear Wang et 



al. (2006)] — the shear banding seen here is thus simply due to the top plate Janiaud et al. 



(2006); Schall and van Hecke (2010)]. We check rheometrically that the top plate drag is given 
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Figure 4: (a) Top view of linearly sheared monodisperse foam. Arrows indicate direction of shear 
at the boundaries, (b) Bubble tracks for monodisperse foam in (a): bubbles slide along planes thus 
move affinely; the resulting velocity profiles are strongly localized (red curve), (c) Bidisperse foam, 
and its corresponding bubble tracks (d): swirly, non-affine motion is clearly visible, and the velocity 
profile exhibits less localization. Ingredients into the model describing these flows: (e) Viscous 
friction between to ordered bubble rows (see inset) as a function of Au. Solid line ~ A?;^/^. (f) 
Stress-strain rate relation for disordered two-dimensional foam: black solid line indicates fit to the 
Herschel-Bulkley expression: a = ay + ^7^, with (5 = 0.36. This power law scaling is highlighted 
by subtracting the yield stress ay (blue data). 



by the classic result of Bretherton F^y^jj 



,2/3 



Bretherton (1961)]. We then extract the scaling of 



interbubble friction with velocity gradient by fitting the model to the velocity profiles. 

Our model predicts rate independent velocity profiles, as seen in the ordered foam, only if the 
interbubble viscous friction scales in the same way as the bubble- wall friction: F^ ~ At'^/'^, where 
At; is the velocity difference between neighboring bubbles. We verify this scaling with rheological 
measurements on two ordered bubble layers gliding past each other, and find that the bare viscous 
friction between bubbles moving past each other indeed scales in exactly the same way with the 
velocity difference between bubbles as the frictional force between bubble and top plate, see Fig. 4e. 
Hence, the rate-independence of the observed flow proflles for ordered foams indicates that the local 
viscous drag law immediately and trivially sets the global drag forces. 
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For the disordered foams, the remarkable thing is that their rate dependent profiles are consis- 
tent with the force balance model if the exponents for bubble-bubble and bubble-wall friction are 
different. In fact, we excellently fit all our velocity profiles by our model with a bubble- wall friction 
as dictated by Bretherton's result and an average interbubble friction that scales as Av^'^^. Since 
the local friction law between bubbles is independent of the bubble configuration, we conclude that 
in this case the translation from local to global drag forces is highly nontrivial. 

We have further checked this remarkable result by bulk rheometry on disordered two dimensional 
foam layers. Bulk foam fiow curves are commonly fit with the phenomenological Herschel-Bulkley 
constitutive relation, linking the stress a with the strain rate 7 : 

a = ay + Aj^ (1) 

with ay the yield stress. In this measurement (Fig. 4f), the flow exponent /3 = 0.36, consistent with 
our indirect determination of /3 via the flow profiles. We thus find that disorder changes the effective 
viscous friction between bubbles in a highly nontrivial way: the averaged viscous dissipation inside 
the foam is enhanced with respect to the ordered case, leading to less localized and rate dependent 
fiows. 

In the next section, we will present a theoretical model that fully captures the nontrivial trans- 
lation from the local exponent of 2/3 to the global exponent of 0.36. 

We also point out here another difference between the measured rheology of disordered foams 
and ordered foams: the disordered foams have a finite yield stress (cy), while the ordered foams 
exhibit pure power law scaling. As we will discuss in section V.C, this is deeply connected to the 
relation between fluctuations and energy dissipation in the system. 

Fluctuations — We have also studied the bubble displacements in detail via tracking of the 



irregular bubble motion [Mobius et al. (2010)]. After subtracting their average velocity, we have 
probed their remaining erratic motion. Since we can probe local strain rates that span multiple 
decades, we can plot the mean squared displacement (MSD) of the bubbles as a function of local 
strain rate. For all strain rates we see that the MSB's cross over from superdiffusive to diffusive 
behavior at a well defined relaxation time U that corresponds to a MSD of (0.14((i))^, which is 



remarkably similar to the Lindemann criterion for cage breaking in colloidal suspensions iBesseling 



et al. (2007)] 



When we extract this for different local strain rates, we observe that it does not scale with 
the inverse local strain rate 7~^, which could have expected expected since both rearrangements 
and fluctuations are entirely shear induced. In contrast, we find that tr ~ ^-o.66±o.05_ ^j^j^^g result 
implies that fluctuations increase for slower flows, in contrast to the commonly accepted viewpoint 
that foam flows tend to quasistaticity when lowering the strain rate, reflected in the expected scaling 
tj. ~ 7~^. This result can be made more explicit by replotting the MSD curves as function of the 
accumulated strain — the lower the strain rate, the larger the MSD at given strain [Mobius et al. 



(2010)]. Hence fluctuations become stronger the slower the flow. 



5 Disorder and fluctuations in foam flows: Simulations and Theory 

In this section we will present a jamming-inspired, theoretical perspective on the flow of foams. 
For completeness, we first introduce several variants of the "bubble model" , a microscopic model 
suitable to simulate the flow of wet foams — readers familiar with, or not interested in the details of 
these models can skip this section. We then discuss the main rheological features that these models 
display in direct numerical simulations. We then show that a basic energy balance equation implies 
that the relative strength of fluctuations grows for decreasing flow rates, and actually diverges in 
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the limit of zero flow rate - consistent with the experimentally observed growth of the diffusivity 
discussed above. We finally outline the main contours of a model for the rheology of soft materials, 
such as foams, near jamming. This model combines recent insights on elasticity near jamming with 
our observations on the nature of the fluctuations, and predicts several scaling regimes that should 
be, and partly have been, observed in the flow of foams. 

5.1 Computer Models for Foam Flow 

There are various versions of microscopic models suitable for simulating the flow of wet foams (and 
other disordered media) . All version we will discuss here stem from the "bubble model" introduced 
in 1995 by Durian [Durian (1995)]. In these models, all bubbles are represented by soft spherical 



(disk-like in 2D) particles, that only interact when in contact. Following Durian, it is typical to 



et al. I p999) ; ,01sson and Teitel 



(2007 


); 


Tighe et al. ( 


2010) 



an appropriate approximation for 
slow flows in which damping dominates inertia. This formally corresponds to the limit in which 
bubble masses are set to zero, and hence at any instant the elastic and viscous forces on each bubble 



balance. Other authors elect to retain inertial terms in the equations of motion iHatano (2008); 



Langlois et al. (2008); Otsuki and Hayakawa (2009)]; in the underdamped limit this produces a 



qualitatively different rheology. 




Figure 5: (a) Elastic and (b) viscous forces between bubbles in the bubble model: when bubbles 
overlap, the strength of their mutual repulsion is a function of the overlap — the bubbles act as 
one-sided springs. The viscous forces are taken to be a function of the velocity difference between 
bubbles sliding past each other. 



To specify a particular variant of the bubble model one must describe the elastic and viscous 
forces between bubbles. Elastic forces are typically taken to be generated by harmonic "one-sided 
springs", see Fig. 5a. These springs have rest lengths equal to the sum of the contacting bubbles' 
radii; their one-sidedness refers to the fact that they apply a force only when the spring is com- 
pressed, thereby guaranteeing that the elastic forces are purely repulsive. Another common force 



law is the Hertzian interaction, which describes emulsions rather than foams [Lacasse et al. (1996)]. 



A "Hertzian spring" is also one sided, but its repulsive force grows as the 3/2 power of the spring's 
compression. In general, then, 

/el ~ <5°=' , (2) 

where 6 is the dimensionless compression of the spring and Ogi = 1 (3/2) corresponds to harmonic 
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(Hertzian) forces — for the simulations we describe below a^i = 1- 

Similarly, it is natural to invoke a viscous force law in which the relative velocity of contacting 
bubbles is damped (Fig 5b): 



f- 

J VI 



vise 

vise 



6ll(A?;ll 
b^{Av- 



(3) 



For the simulations we describe below Ovisc = Ij which corresponds to linear viscous dissipation. 
Real foams, however, are believed to have Ovise < 1, with the actual value depending on the 

D09 )] . For the mobile surfactants used in the experiments described 
our measurements indicate avisc = 2/3, see Fig. 4(e). We shall see that changes in the 



surfactant used |Denkov et al. 
above, 

microscopic exponent Ovise have measureable ramifications for the macroscopic rheology. We also 
take 6" = 6"*", so that the dissipation is indifferent to whether the particles are moving together, 
moving apart, or sliding. Other authors have taken 6-*- = 0, so that sliding does not dissipate 
[Hatano | (2010); Otsuki and Hayakawa (2009)], though this case may be pathological. 



Finally we point out that there is another, qualitatively different, viscous force law commonly 
used in the literature iDurian 



(1995 


1997 


); 


Ono et al. 


( 


2002 


); 


Tewari et al. 


( 


1999 


); 


Olsson 



and Teitel (2007)]. Known as "mean field" dissipation, it specifies that each bubble experiences 
a damping force proportional to the difference between its instantaneous velocity and the mean 
velocity t/iin = 7X y of a bubble at its position: 



fj 



MF 



-buFiv-Vlir 



(4) 



This expression assumes the linear velocity profile established by Lees-Edwards boundary conditions. 
It sacrifices more realistic modeling for numerical convenience: in overdamped dynamics, Eq. ([4]) 
can be simulated far more efficiently than Eq. ([3]). 



5.2 Phenomenology 

In the static limit, all microscopic models introduced above reduce to the simple models that are 
studied to probe elastic properties near the jamming transition. Most crucial, in order for the 
material to have a finite rigidity, they need to be packed at a sufficiently large packing fraction 
— the critical packing fraction corresponds to extremely wet foams where the bubbles just touch. 
Once such packings are created, they are subjected to a fixed strain rate and the shear stress can 
be measured numerically — in some cases, the stress has been fixed and the strain rate has been 
measured 



Olsson and Teitel 



(2007)], but in all cases, flow profiles have been found to be simply 
linear (i.e., no shear banding — when we include additional drag terms, modeling interactions with 
a top plate, the same shear banding that is seen experimentally is recovered). 

There are certain generic features of the macroscopic rheology that prevail for </> > in sim- 
ulations of any version of the bubble model discussed above. Namely, the constitutive relation is 
qualitatively similar to the Herschel-Bulkley constitutive relation defined in Eq. ([!]) The exponent 
f3 is typically less than one, a property known as "shear thinning." In numerics one observes that 
ay varies with (f) and vanishes at (^c, which motivates us to write ay ~ Acp^ for some A. Note 
that because a = ay in the limit 7 — )■ 0, it is natural to expect that the exponent A depends on 
ttei but not on avisei a dynamical quantity. By contrast /3 is a dynamical exponent which will turn 
out to depend on both Ovise and a^- The challenge is to identify how A, A, and /3 depend on the 
microscopic details of the model. 

^We note here that we do not beUeve that the flow curves are exactly of this form, although this is almost always 
a good approximation 
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5.3 Fluctuations 



Slow foam flows are dominated by fluctuations. For example it is widely believed that the main 
contribution to dissipation in dry foams stems from localized neighbor-swapping, or Tl, events 



[Hohler and Cohen- Addad (2005)]. Here we focus on the wet limit, where Tl events are not well 
defined. Nevertheless, here again there is a precise relation between the fluctuations and dissipation. 
Power balance — If one shears a foam (or any other system) at a rate 7, and it resists this flow 
with a shear stress a, the product of the two gives the average rate of work done on the system. This 
work has to be dissipated, and for wet foams the only dissipative mechanism is the relative motion 
of neighboring bubbles. Let us first take a simple model, where the viscous force fvisc between 
bubbles is proportional to their relative velocity Av. In a time interval At, the amount of energy 
dissipated per pair of sliding bubbles is then fvisc x At x Av, which scales as Av'^At. The crucial 
observation is that the energy fed into the system by shear needs to be balanced by the amount of 
energy dissipated by local bubble sliding: 

^7 ~ {Av^) , (5) 

where the brackets denote averaging over the system]^ 

We now show the powerful consequences of this simple relation. Let us first imagine that our 
foam behaves as a Newtonian fluid, for which a = rj'j. In that case the steady state power balance 
equation Eq. ^ requires that the rj-j"^ ~ (^Av'^'j, so that \Av\ and 7 scale in the same way. 

However, foams do not behave as Newtonian fluids, particularly for slow flows. Let us assume 
that the stress varies as a power law of the strain rate, which is a good approximation of a typical 
Herschel-Bulkley rheology at high flow rate: cr ~ 7''. Substituting this into the steady state 
power balance equation, Eq. ([5|, one obtains 7-^+^ ~ (^At;^^, so that \Av\ and 7 do not scale in 
the same way; rather lA^I ~ j^^^^y^. For the typical case that /3 < 1, one concludes that the 
velocity fluctuations decay sublinearly with strain rate, or in other words, that the relative velocity 
fluctuations, |At;|/7 diverge as 7(^~^)/2. This divergence is even stronger for slow flows for which 
the stress reaches a plateau value: a ^ ay. The steady state power balance equation Eq. ^ then 
predicts IAvI/t- to diverge as 7"^/^. 

The divergence of the fluctuations suggests that there is no well-deflned quasistatic limit, at 
least not as far as the trajectories of the bubbles are concerned: the slower you go, the bigger the 
fluctuations. Imagine you are given two movies of foam flows at two different flow rates — but 
you are not told which one is the fastest, nor do you know the frame rate of these movies. By 
adjusting their playback speed, you can make the average flow rate of the two movies equal — the 
energy balance argument predicts that the amount of fluctuations would be largest in the movie 
corresponding to the slowest flow. We stress here that, at least qualitatively, this is precisely what 
we observed in the experiments on diffusion of foam bubbles in flow. 

Characteristic Scales — The fluctuations introduce a nontrivial dynamical time scale into 
the problem: tdyn = f^/l^^^l) where d is a typical bubble diameter — this time can be thought of as 
characterizing the time scale over which local rearrangements take place. In the simplest case of a 
Newtonian fluid, for which |At;| ~ 7, tdyn is nothing more than the inverse strain rate. 

For non-Newtonian fluids, however. At; scales nontrivially, hence the dynamical timescale does, 
as well. For a pure power law fluid we find tdyn ~ 7(1+/^)/^^ and for the yield stress case we find 

One can translate this time scale into a characteristic strain scale necessary to induce rearrange- 
ments: 7dyn = ^dyn7- This characteristic strain thus vanishes for slow fiows — the slower the fiow, 

^In all scaling arguments that will follow we focus on the typical scale of quantities, and ignore correlations — we 
do not believe these will change our results in an essential manner. 
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the smaller the overall strains necessary to induce substantial rearrangements. 

Relation between local viscous drag, fluctuations and rheology — Let us now clarify 
the relations between local dissipation, fluctuations and global rheology. First, let us generalize the 
viscous force to be nonlinear: /vise oc |Az;|". Now assume that the macroscopic rheology is of the 
form (T oc 7^^. The power balance equation then reads 7'^+-'^ ~ (^Af^"^"), which will give nontrivial 
fluctuations whenever a 7^ /3. 

This shows that details matter in the sense that different local dissipative laws directly affect 
the precise form of the balance Eq. ([s]) — in this sense, the physics is not universal. On the other 
hand, as long as /3 < a, the fluctuations become stronger for slower flow. We will see that this 
inequality is in general satisfied, so that the divergence of fluctuations is robust. 

The balance arguments can also be used to rationalize part of our experimental findings. Recall 
that in rheological experiments in which two ordered rows of bubbles were slid past each other with 
velocity Av Katgert et al. (2008, 2009)], see Fig. 4e, the time-averaged shear stress a was found 
to scale as a simple powerlaw: a ~ Af". In particular, there was no force plateau at low Av. In 
contrast, for a disordered bubble raft, the stress had a Hershel-Bulkley form with a plateau for low 
strain rates and a power law different from a: o" ~ ctq + Aj^, see Fig. 4f. 

The finite plateau we can understand as follows. As we have seen, the energy balance argument 
shows that for flows in the regime where the stress is on a plateau, the relative fluctuations must 
diverge when the strain rate goes to zero. Since the fluctuations are constrained in the ordered 
system, they cannot diverge, energy cannot be dissipated strongly enough, and as a result there 
cannot be a yield stress — consistent with our findings. In the disordered system, the fluctuations 
are not constrained, and nothing forbids the emergence of a flnite yield stress. 

We stress here that the role of disorder is to facilitate the nontrivial role of the fluctuations. In 
ordered systems, such as the ordered foams we discussed above, whole rows of bubbles slide past 
each other, and there is simply not enough freedom to allow for large fluctuations. In disordered 
systems, the bubbles have much more freedom to choose their path. For the case of elasticity 
of disordered media, discussed in section II, disordered bubble motion was intimately connected 
to anomalous scaling of the shear modulus. Here we see that for the flow of disordered media, 
nontrivial fluctuations are connected to anomalous scaling of the stress-strain rate relation. 

The difference between the local and global rheological exponents a and /3 is also intimately 
related to the fact that the fluctuations are rate dependent. But the single balance equation Eq. ^ 
is not sufficient to predict both /3 and the fluctuations — we need additional arguments to obtain a 
set of closed equations. In the next section we will introduce such equations, and produce a definite 
prediction of /3 as function of the local drag force a. 



5.4 Scaling model for foam flows 

We now introduce a model that seeks to explain bubble model rheology for (j) > cpc- Our goal is 
to give the fiavor of the model, without delving too far into details. For simplicity, we first fix the 
elastic (oei) and viscous exponent (avisc) to one. 

The scaling model has three ingredients. These are 

1. the system is in power balance 

2. a flowing foam can be mapped to a static system that has been sheared through an effective 
strain 7efj 

3. the stress as a function of 7efr is given by the constitutive relation for sheared disordered spring 
networks 
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One of the difficulties in explaining this model is that there are several different regimes with 
different forms for the effective strain and stress-strain relation. We first will describe the model in 
two of the simplest limits, and then briefly point out how all these results can be generalized. 

Critical Regime — In the limit of very wet foams, i.e., (/> — >• (pc the model is particularly 
simple. We postulate that the effective strain in the system then simply equals the dynamical 
strain introduced above: 7dyn ~ 7/|Af |. For the equation for the stress, we build on recent work 
that shows that near jamming, the linear shear modulus vanishes and the stress scales as cr ~ |7|7- 
[Wyart et alT~| ( |2008D ]. Now we have a closed set of three equations (we drop all absolute values 



and averages): 

0-7 ~ Av^ . (6) 

7cfr ~ 7/ Aw . (7) 

a ~ 7,2^. (8) 

We have identifled 7 with 7efr; this is point 3 in the above list. These equations can easily be solved 
and yield a prediction for the rheology of the form a ~ j^^"^, which is in very good agreement with 



numerical simulations of the bubble model [Tighe et al. (2010)] 



Yield Stress Regime — In the limit of very slow flows and for (j) > the model also becomes 
very simple. In that regime, the effective strain is no longer expected to be dominated by dynamic 



effects, and we postulate that the effective strain in the system then equals A0 := — 0c Tighe et 



al. (2010)]. For the equation for the stress, recent work shows that away from jamming, the shear 



modulus G scales with the distance to jamming as G ~ Az ~ yA^- In summary: 

(T7 ~ Aw^ . (9) 
7efr ~ A0. (10) 
a ~ A/A07eff . (11) 

These equations can easily be solved, and yield as prediction that the stress is a constant oc 

Transition Regime — The general form of the model is obtained by combining these regimes, 
so that 7eff = aiA(/) + a27/Af , and a = a■i^J~K4>le^^ + o^47efr) where Oj are numerical constants to be 
determined. It turns out that there are three different regimes, and in all cases, either the flrst or 
the second term in these sums dominates |^ The third regime we refer to as the transition regime, 
and here the three equations are: 

0-7 ~ At;^ . (12) 
7efr ~ 7/At; . (13) 
a ~ \/A07eff , (14) 

which yields that the stress scales as cr ~ (A(/))^/^7^/^ — this third regime came as a complete 
surprise, and its existence implies deviation from the usual Herschel-Bulkley phenomenology. This 
deviation is difficult to observe in numerics but consistent with our observation that very often, 
Herschel-Bulkley-flts to experimental or numerical data underestimate the data in the crossover 
regime between yield stress and power law behavior - precisely the regime where the transition 
rheology is predicted. 

By substituting the various solutions we have obtained into the general expressions for the 
effective strain and stress, one can do a self consistency check to see for which strain rates, which 

^In principle there are four regimes, but the combination of terms in this fourth regime never dominate the physics 



see Tighe et al. (20101] 
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regime should dominate — we find that the Yield Stress regime dominates for 7 < IS.i^l'^ ^ the 
Critical Regime for 7 > A(/>^ and the Transition regime in the range in between — which can span 
arbitrarily many decades when A(/) tends to zero. 

General Microscopies — What happens when the microscopic exponents (ad) and (ovisc) are 
unequal to one? 

For the Yield Stress regime, the stress-strain relation is affected by the value of Oeh as it sets 
the scaling of the shear modulus G with Ac/). Prior work has shown that for nonlinear interactions, 
the stress strain relation becomes a ~ A(/»"<=' — 1/2 [O'Hern et al. (2003)]. Thus the yield stress 
varies with Uei according to A = Ogi + 1/2- For Oei = 3/2, believed to describe dense emulsions and 

(|2010a[)], this gives ay ~ Acj)^ 



Lacasse et al. 



(1996); 



Nordstrom et al. 



(2010) 



microgel suspensions 
good agreement with experimental measurements by Mason et al. 
which both find A 



m 



(1996) and Nordstrom et al. 



2. 



For the Critical Regime, both the elastic and viscous exponent matter. Let us keep the elastic 
exponent equal to one, which is a good approximation for foams. Recall that f3 should, in general, 
depend on Ovisc- Repeating the analysis in the critical rGginic for arbitrary o^visc 

yields 



/3 



2avisc 

Ovisc + 3 



(15) 



For Ovisc = 2/3, as is the case for mobile surfactants, we obtain (3 = 4/11 ~ 0.36, in remarkably 
good agreement with the experimentally determined value /3 ~ 0.36 inlKatgert et al. 



(2008). 



Finally, we point out here that for dense granular media, the main dissipation comes from sliding 
friction, which could be seen as a viscous interaction with exponent zero — the frictional forces do 
not change appreciably with the sliding rate. In granular media, it has been taken as a triviality that 
both the local and global interactions are very similar, namely frictional. We doubt that granular 
flows are affine, but we believe that this correspondence between local and global flow behavior is 
a lucky coincidence - we note here that for ctyisc = 0, the global flow exponent /3 also becomes zero 
in our model. 



5.5 Critical collapse 

A powerful test of any critical scaling prediction is to plot numerical or experimental data in rescaled 
coordinates and look for collapse to a master curve. The quality of collapse is a strong test of theory, 
and the master curve itself is a revealing depiction of the underlying physics. The rheological model 
described above indeed predicts such master curves, but there is some subtlety owing to the existence 
of not one but two crossover strain rates. 

To illustrate the subtlety, let us first consider an example with only one crossover. The Herschel- 
Bulkley flow rule contains two regimes, a yield stress regime a ~ ciy ~ Acj)^ and a critical regime 
£7 ~ 7^. The crossover between the two regimes occurs at the strain rate 7* for which the two 
regimes are comparable: 7* ~ A(/)^/^. Dividing by A(p^, the Herschel-Bulkley flow rule can be 
rewritten ^ 

^ const + a( \ ] (16) 



A(j)^ \A^^/I^ 



provided A does not depend on Acp as well Katgert et al. ( 2009 )] . Note that the term in parentheses 
is proportional to 7/7*, and that all the dependence on Acj) is contained in the ratios a/Acj)^ and 
7/A0^/^. Therefore if the rheology is described by a Herschel-Bulkley rule (or another flow rule 
with a single crossover rate), plotting the rescaled stress and strain rate, 

^= A^ ^"*^^= A^b' ^1^) 
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Figure 6: Scaling collapse of flow curves obtained from bubble model simulations as described in 
Section 5.1. Flow curves are obtained for 4 decades in strain rate and 3 decades in A, see legend. 
Rescaled coordinates are a j and 7/A0^, appropriate for parameters spanning the Transition 
and Critical regimes. Dashed lines are guides to the eye with slopes 1/3 and 1/2. Inset: Boundaries 
between the yield stress, transition, and critical regimes in the — 7 plane. 



should produce collapse to a master curve of the form const + Ax^ . Alternatively, if the exponents 
A and /3 are unknown, such a scaling plot can be a way to determine them: one looks for the values 
of A, /3, and that produce the best data collapse. 

Such scaling plots can be made for any curve characterized by two qualitatively different regimes: 
the first step is always to identify the crossover strain rate (or whatever quantity is on the x- 
coordinate). If a curve has more than two regimes, and hence more than one crossover, one generally 
cannot collapse all data by resaling with the single scaling parameter Ai;^. The flow model from 
above predicts not two but three scaling regimes, hence, according to the model, one should not 
expect to be able to make a scaling plot that collapses all data near for a broad range of 7. 
One can, however, still collapse data from any two adjacent regimes. Consider the Transition and 
Critical regimes for ad = 1 and Ovisc = 1, for which a ~ ^^I'^^^l'^ and a ~ ^'^^'^ ^ respectively. 
The crossover rate is 7* ~ Ac/)^, and data from these two regimes should collapse for the rescaled 
coordinates 

o- , ■ 7 

— - and 7 = — — „ 

This is shown in Fig. [6] and the collapse is indeed very good. Similarly, data from the yield stress 
and transition regimes can be collapsed by plotting 

^=Aiv^ ""^^=A^- ^''^ 

Note, however, that generating data in the yield stress regime is numerically challenging; Fig. [6j 
for example, does not have enough data points in this regime to allow for a convincing test of the 
model. 



~^=—^. and ^ = — 2 . (18) 
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Figure 7: Numerical solution to our scaling model (♦) and Herschel-Bulkley model with the same 
asymptotic behaviour (dashed line) compared. The Herschel-Bulkley model underestimates the 
cross-over as it ignores the Transitional regime. 



It is interesting to ask what happens if we attempt to make a scaling plot that assumes a 
single crossover, as in the Herschel-Bulkley flow rule, when in fact there are multiple crossovers. In 
effect, this ignores the transition regime, and one should find that data from the transition regime 
fails to collapse. Fig. [7] shows a plot of data generated by numerically solving our scaling model 
(aei = Ovisc = !)• The dashed line indicates a Herschel-Bulkley law: ay + Aj^^'^, with ay and A 
chosen so as to match the asymptotics of the model. Note that the curve has a soft elbow: the 
crossover is much slower than the Herschel-Bulkley form would predict. While real experimental 
data likely are noisy enough to mask the poor collapse in this regime, the soft elbow does appear 
to be a feature of real flow curves: fits to the HB flow rule typically fare poorly near the crossover 
point and underestimate the stress, consistent with the existence of a transition regime. 



5.6 Conclusion and Outlook 



In this paper, we have proposed that static and dynamic properties of foams can be captured by 
recent ideas stemming from jamming - and we have also illustrated how some of these theoreti- 
cal ideas grew out of experimental questions ("how does a foam lose rigidity") and observations 
(" fluctuations in foams seem to grow" ) . 

We suggest two broad directions for future research. On the one hand, there is a general lack 
in experimental data confronting some of the predictions, for example in the degree of nonaffinity 
as a function of wetnes s [van Hecke (2010)], or the distribution of relative velocity fluctuations as 
function of strain rate Tighe et al. (2010)]; more experiments are called for. In particular there 
are only a handful of studies of the mechanics and flow of very wet foams or emulsions - the recent 



advent of combined 3D confocal imaging and rheology of emulsions [Zhou et al. (2006); Clusel 



et al. (2009); Paredes et al. 



(2011)] should open new experimental avenues to probe the critical 
regime. 

On the other hand, many of the commonly observed phenomena in soft materials are without a 
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sound description. Here we have focused on steady state rheofogy, but in practice, oscillatory rheol- 
ogy is the preferred experimental tool to capture the general visco-elastic behavior of soft materials. 
One of us, has recently developed a theoretical description of linear oscillatory rheology [Tighe 



(2011)], but there are many other open questions. What is the physics of plastic rearrangements? 
What is the fate of Tl evens when the foam gets increasingly wet? What about memory effects and 
reversibility [Lundberg et al. (2008)]? 



We look forward to many more surprises in the rich physics of collections of bubbles and droplets. 
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